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Abstract 



We present the non-Abelian Stokes theorem for the Wilson loop in various forms and 
. discuss its meaning. Its validity has been recently questioned by Faber, Ivanov, Troit- 

skaya and Zach. We demonstrate that all points of their criticism are based on mistakes 
in mathematics. Finally, we derive a variant of our formula for the Wilson loop in lat- 
tice regularization. 



1 Introduction 

One of the main objects in the Yang-Mills theory is the Wilson loop or holonomy; it is 
defined as a path-ordered exponent, 

W r = —- Tr P exp i idr — A a u T a , (1) 
d(r) J dr ^ 

where x M (r) with < r < 1 parametrizes the closed contour, A ^ is the Yang-Mills field (or 
connection) and T a are the generators of the gauge group in a given representation r whose 
dimension is d{r). 

It is generally believed that in three and four dimensions the average of the Wilson loop 
in a pure Yang-Mills quantum theory exhibits an area behavior for large and simple contours 
(like flat rectangular). This should be true not for all representations but those with W-ality' 
nonequal zero; in the simplest case of the 577(2) gauge group these are representations with 
half-integer spin J. 



One of the difficulties in proving the asymptotic area law for the Wilson loop in half- 
integer representations (and proving that in integer representations it is absent) is that the 
Wilson loop is a complicated object by itself: it is impossible to calculate it analytically 
in a general non-Abelian background field. Meanwhile, it is sometimes easier to average a 
quantity over an ensemble than to calculate it for a specific representative. However, in case 
of the Wilson loop the path-ordering is a serious obstacle on that way. 

A decade ago we have suggested a formula for the Wilson loop in a given background 
belonging to any gauge group and any representation [|TJ . In this formula the path ordering 
along the loop is removed, but at the price of an additional functional integration over all 
gauge transformations of the given non-Abelian background field. This formula is discussed 
below, in section 2. Furthermore, the Wilson loop can be presented in a form of a surface 
integral ||, see section 3. We call this representation the non-Abelian Stokes theorem. It is 
quite different from previous interesting statements |3|, [|, |5], || also called by their authors 
'non-Abelian Stokes theorems' but which involve surface ordering. Our formula has no 
surface ordering. A classification of 'non-Abelian Stokes theorems' for arbitrary groups and 
their representations has been given recently by Kondo et al. who used the naturally 
arising techniques of flag manifolds. 

Though these formulae usually do not facilitate finding Wilson loops in particular back- 
grounds, they can be used to average over ensembles of Yang-Mills configurations, and in 
more general settings, see e.g. |8|, |, 0, [UJ. 

Our formula for the Wilson loop have been recently questioned by Faber, Ivanov, Troit- 
skaya and Zach (FITZ) [|ll|]. In section 4 we show that all points of their criticism are due to 
mistakes in mathematics, which we thoroughly locate, one by one. An alternative formula 
for the Wilson loop proposed by FITZ is also mathematically inconsistent, and we pinpoint 
their concrete errors. 

Finally, in section 5 we present a variant of our formula for the Wilson loop in lattice 
regularization. It appears to be very similar to that presented in section 2. 



2 Formula for the Wilson loop 

Let r parametrize the loop defined by the trajectory x^{r) and A(t) be the tangent compo- 
nent of the Yang-Mills field along the loop in the fundamental representation of the gauge 
group, A(t) = A^dx^ /dr, Tr(t a t b ) = \5 ab . The gauge transformation of A(r) is 

A(r) -> S-\t)A(t)S(t) + iS-\r)^S(r). (2) 

Let Hi be the generators from the Cartan subalgebra (i = 1, r; r is the rank of the gauge 
group) and the r-dimensional vector m be the highest weight of the representation r in which 
the Wilson loop is considered. The formula for the Wilson loop derived in ref. |I[ is a path 
integral over all gauge transformations S(t) which should be periodic along the contour: 

W r = J DS(t) exp ijdrTr [m t Hi (S^AS + iS^S)] . (3) 

Let us stress that eq. (|3]) is manifestly gauge invariant, as is the Wilson loop itself. For 
example, in the simple case of the SU (2) group eq. (||) reads: 
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Wj = J DS(t) exp i J J dr Tr [r 3 (SU5 + iS ] ' S)] (4) 

where r 3 is the third Pauli matrix and J = |, 1, |, ... is the 'spin' of the representation of 
the Wilson loop considered. 

The path integrals over all gauge rotations (|3],[|D are not of the Feynman type: they do 
not contain terms quadratic in the derivatives in r. Therefore, a certain regularization is 
understood in these equations, ensuring that S(r) is sufficiently smooth. For example, one 
can introduce quadratic terms in the angular velocities iS*S with small coefficients eventually 
put to zero; see ref.|lj] for details. In ref. U eq. (£|) has been derived in two independent ways: 
i) by direct discretization and ii) by using the standard Feynman representation of path 
integrals as a sum over all intermediate states, in this case that of an axial top supplemented 
by a 'Wess-Zumino' type of the action. Another discretization but leading to the same result 
has been used recently by Kondo j?| . A similar formula has been used by Alekseev, Faddeev 
and Shatashvili |L2| who derived a formula for group characters to which the Wilson loop 
is reduced in case of a constant A field actually considered in [T2[]. In ref.[|I3[] eq. (]3p has 
been rederived in an independent way specifically for the fundamental representation of the 
SU (N) gauge group. Finally, at the end of this paper we consider the Wilson loop in lattice 
regularization and derive a formula similar to eq. (|j). 

The second term in the exponent of eqs. (|3],|4]) is in fact a Wess-Zumino '-type action, 
and it can be rewritten not as a line but as a surface integral inside a closed contour. Let 
us consider for simplicity the SU(2) gauge group and parametrize the SU(2) matrix S from 
eq. (|j) by Euler's angles, 

S = exp(— ryr 3 /2) exp(-i/3r 2 /2) exp(-iar 3 /2) 



,- °±3 • j a ~"f 

cos ~ e 2 — sm | e 2 



P J 



sin^e 2 cos | e 2 



(5) 



The derivation of eq. @ implies that S(t) is a periodic matrix. It means that a ± 7 and (3 
are periodic functions of r, modulo Atc. 

The second term in the exponent of eq. (|]) which we denote by $ is then 

$ = J drTr(T 3 iS^S) = J dr d(cos j3 + 7) 

= J dr [d(cos/3 — 1) + (a + 7)] = J dr a (cos (3 — 1). (6) 

The last term is a full derivative and can be actually dropped because a + 7 is 47r-periodic, 
therefore even for half- integer representations J it does not contribute to eq. Notice 
that a can be 27r-periodic if 7 (which drops from eq. ©) is 2tt, 6tt, . . .-periodic. If cn(l) = 
a(0) + 2ir k we shall say that «(t) makes k windings. Integration over all possible a(r) 
implied in eq. (§) can be divided into distinct sectors with different winding numbers k. 

Let us prove that eq. (|6|) can be written as an integral over any surface spanned on the 
contour (we shall call it a 'disk'), 

$ = - J drda e abc Uj n a ftn fe ^n c , i, j = r, a, (7) 
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where n is a unit 3-vector, 



n a = -Tr (S^r a Sr 3 ) = (sin (3 cos a, sin (3 sin a, cos/5). (8) 

It is understood that n or a and /3 are continued inside the disk. We denote the second 
coordinate by a such that a = 1 corresponds to the edge of the disk coinciding with the 
contour and a = corresponds to the center of the disk. Let us establish what conditions 
should the continuation n(cr, r) satisfy. 
Eq. (|7|) can be identically rewritten as 

$ = J J drda [d T ad a (cos (3 — 1) — d a ad T {cos (3 — 1)] 
= J J drda {d a \d T a(cos(3 — 1)] — d T [d a a(cos/3 — 1)]} 
= / dr [d T a (cos (3 - l)] a Zl - f da [d a a (cos (3 - 1)]^ . (9) 

The second term here is zero, for the following reasons. Let ot(r) belong to the sector 
with winding number k. The periodicity property of a(a, r) cannot change abruptly as we 
continue it inside the disk. We can write a(a, r) = 2irkT + a(a, r) where a is strictly periodic, 
as well as cos (3. Therefore, the integrand of the second term in eq. @ is zero by periodicity. 

Let us now consider the first term in eq. @. We want the surface integral (0) to reproduce 
the contour integral (((]), up to a possible contribution of 4irl with integer I, coming from the 
center of the disk, o = 0. Such a contribution is allowed since it does not affect eq. (f|) even 
for half-integer J. The contribution of the first term in eq. @ at a — is 

- / dT[2Trk + d T a(0,T)][cosl3(0,T) - 1] = 2?rA;[cos/?(0) - 1] =4ttZ, (10) 
Jo 

where we have used that a is periodic while (3(0, t) is in fact independent of r, otherwise 
n 3 = cos (3 would be not uniquely defined at the center of the disk. Eq. fllCf ) means that 
either a(r) belongs to the winding number k = sector, then (3(0) is arbitrary, or, if k ^ 
then (3 = 0, it meaning that n 3 = ±1 at the center of the disk. Notice that for \n 3 \ ^ 1 
the components n 1 ' 2 are nonzero and are varying very rapidly as function of r at the point 
cr = 0. These conditions can be summarized as the condition that n(cr, r) should be smooth 
and uniquely defined inside the disk. 

Let us note that if the surface is closed or infinite the r.h.s. of eq. (^) is the integer 
topological charge of the n field on the surface: 



Q = _L / dadT e abc eij n a din b d jn c . (11) 



1 

8^ 

Eq. (0) can be also rewritten in a form which is invariant under the reparametrizations 
of the surface. Introducing the invariant element of a surface, 

d^ = dadr {^4- - d 4- d -^\ = <T rf(Area), (12) 



dr da dr da 



one can rewrite eq. (0) as 
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$ = I J d 2 S^e abc n a d^n b d u n c . 



We get for the Wilson loop 



W, 



Dm 



a, t) exp 



U J dr(A a n a ) + y J d 2 S^e abc n a d ll n b d u n c 



(13) 



(14) 



The interpretation of this formula is obvious: the unit vector n plays the role of the 
instant direction of the color 'spin' in color space; however, multiplying its length by J does 
not yet guarantee that we deal with a true quantum state from a representation labelled by 
J - that is achieved only by introducing the 'Wess-Zumino' term in eq. ([14]) : it fixes the 
representation to which the probe quark of the Wilson loop belongs to be exactly J. 



3 Non-Abelian Stokes theorem 

We can now rewrite the exponent in eq. (1^) so that both terms appear to be surface integrals 
& 

Wj = J Dn(a, r) exp y J d 2 S^ (-F> a + e abc n a (D,n) b (D u n) c ) , (15) 
where D ab = d^8 ah + e acb 'A^ is the covariant derivative and F* =d^A a v - d v A a u + e abc A\ A c h 



is the field strength. Indeed, expanding the exponent of eq. (15) in powers of A^ one observes 
that the quadratic term cancels out while the linear one is a full derivative reproducing the 
A a n a term in eq. (|14"D; the zero-order term is the Wess-Zumino' term ([?D or Note that 
both terms in eq. (|l^) are explicitly gauge invariant. We call eq. fli~5|) the non-Abelian Stokes 
theorem. We stress that it is different from previously suggested Stokes-like representations 
of the Wilson loop, based on ordering of elementary surfaces inside the loop || |], ||, ||. It is 
understood that the functional integration measure in eq. flI5D is such that Wj = 1 for zero 
field. 

The gauge-invariant field strength appearing in eq. (|T§), 

G, v = F« u n a - e abc n a (D^n) b {D v n) c , (16) 



coincides in form with the field strength introduced by Polyakov ]14| and 't Hooft |L5] in 
connection with monopoles. In that case the unit-vector field n a has the meaning of the 
direction of the elementary Higgs field, <fr a /\(j)\. 
One can introduce a 'monopole current', 



3*u — Tr e ^un\9uG K \. (17) 

P 87T 



Using equations 



e» VKX DtF b RX = 0, e^ x e abc (D u n) a (D K n) b (D x n) c = 0, 
e^x{D u D K n) = 0, n a (D K n) a = 0, 
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one can easily check that this current is zero. In mathematical language it is the statement 
that any exact two-form is complete. There is a subtlety here, however. Let us consider a 
configuration which we shall call the 'Wu-Yang monopole', n a = x a /r (not to be confused 
with the monopoles of the field). Let us take a contour of unit radius about the origin 
lying in the z = plane, the equator. In terms of angles the monopole corresponds to 
a(r) = = 27tt, (3 = 9, see eq. (§). At the equator 9 = ir/2, and we have from the line form 

$= f dr d T a (cos (3 - 1) = C dr 2n (-1) = -2tt. (18) 
Jo Jo 

Let us now compute $ in the surface form: first, over the upper hemisphere whose edge is 

the equator, second, over the lower hemisphere. For the upper hemisphere we introduce the 

variable a = 1 — cos 8, so that a = at the north pole and a = 1 on the equator. Then 

cos (3 — 1 = — a, and eq. (|7|) reads 

$= f f 1 drda d T (2-KT)d a {-a) = -2vr, (19) 



io Jo 

in agreement with eq. ([18]) . For the lower hemisphere we define the variable a = 1 + cos 9 so 
that cr = at the south pole and a = 1 on the equator. Then cos (3 — 1 = a — 2, and eq. (0) 
reads 

$= f 1 f 1 drda d T (27iT)d a (cr - 2) = +2vr. (20) 
Jo Jo 

There is no arithmetical contradiction here since in deriving the surface form @ we have 
allowed for the contribution from a = to be a multiple of 4n, and it is exactly what has 
happened here. Simultaneously, it is in correspondence with the fact that the integral over 
the full sphere around a monopole (see eq. (O)) is 4ir. In both cases the contribution of the 



'Wu-Yang monopole' configuration to the functional integral ([15]) is 



exp ±2niJ= (— 1) , (21) 

irrespective of whether the surface is drawn above or below the singularity of the n field. 
Notice that eq. ([21]) makes a clear distinction between integer and half-integer representa- 
tions. 

If, however, we take the surface which passes exactly through the center of the monopole, 
we shall be in trouble. For example, choosing the equatorial plane we have n 3 = everywhere 
on the plane, hence the surface integral (0) is zero, in contradiction with eq. (|18D. It is 
because we have drawn the surface directly through the singularity and have thus violated 
the condition needed for the applicability of the surface form (^). Does it mean that singular 
configurations of the type n a = x a /r are altogether forbidden in the functional integral (|T5D? 
Probably not, as long as such singularities have a zero-measure probability to hit the surface. 

Let us integrate the 'monopole charge density' j'q over a volume surrounded by a closed 
surface. According to the Gauss theorem, 

8tt J d%j; = - J d 2 S^e abc n a d^n b d u n c = -8ttQ, (22) 
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where Q is the integer winding number (|Tl"| ) of the n field over the closed 2-surface surround- 
ing the 'monopole'; in the concrete case of the field n a = x a /r one has Q — 1. Since the 
surface can be drawn as close to the center as one wishes, it is tempting to say that the 
'Wu-Yang monopole' has j'q = —5^(r). Since computing j'q involves singular operations, 
one can regularize the singularity, for example by replacing n a — > x a /(r 2 + 5 2 ) 1 / 2 , 5^0. 
This regularization, indeed, leads to jl = -3<5 2 /(r 2 + 5 2 ) 5/2 / (An) — > -5 (3) (r). It should 
be stressed, however, that this regularization violates the condition n 2 = 1 which has been 
crucial to prove that j* = 0. 

To conclude this discussion: The Wilson loop is defined as a contour integral, so when 
one writes it in a surface form one has to take care that it does reproduce the contour 
form. Generally, singular n(a, r) leading, after regularization, to nonzero j* imply that the 
result depends on how one draws the surface, which is unacceptable. However, e.g. a 'gas' of 
singularities of the type n a = x a /r with quantized Air flux (in d = 4 it is a 'gas' of worldlines) 
is still admissible since the probability that the singularity occurs exactly on the surface has 
zero measure. 

Let us briefly discuss gauge groups higher than SU (2): for that purpose we have to return 
to our eq. ([J). Eq. (|3|) is valid for any group and any representation. It is easy to present it 
also in a surface form. We denote the combination of Cartan generators rriiHi where m is 
the highest weight of a given representation r by 7i r . Using the identity 



e {j di Tr HrS^VjS = e tJ Tr H r IS' 1 ( V, VjS) + (S' 1 V,) (VjS) 



€ij Tr H r 



(23) 



we can present eq. (|3]) in a surface form: 

W r = J DS(a,r) exp i J dS^TrH r l-^S^F^S) + (S^V^) (V U S) 



(24) 



Actually, eqs. (0,0) depend not on all parameters of the gauge transformation but only 
on those which do not commute with the Cartan combination 7i r = rriiHi. In the SU(2) 
case one has rriiHi = Jr 3 , J = 1/2, 1, 3/2, . . ., since SU (2) is of rank 1, and there is only one 
Cartan generator. Therefore, in the SU(2) case one integrates over the coset SU(2)/U(1) 
for any representation; this coset can be parametrized by the n field as described above. 

In case of higher groups the particular coset depends on the representation of the Wilson 
loop. For example, in case the Wilson loop is considered in the fundamental representation 
of the SU(N) group the combination rriiHi is proportional to one particular generator of the 
Cartan subalgebra, which commutes with the SU (N — 1) x U(l) subgroup. [In case of SU(3) 
this generator is the Gell-Mann As matrix or a permutation of its elements.] Therefore, the 
appropriate coset for the fundamental representation of the SU(N) group is SU(N) / SU(N— 
1) /U{\) = CP 1 ^^ 1 . A possible parametrization of this coset is given by a complex iV-vector 
u a of unit length, u\u a = 1. To be concrete, the Cartan combination in the fundamental 
representation can be always set to be rriiHi = diag(l,0, ... ,0) by rotating the axes and 
subtracting the unit matrix. In such a basis u a is just the first column of the unitary matrix 
5* while is the first row of S^ . Unitarity of S implies that u^u 01 = 1. 
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In this parametrization eq. (H) can be written as 



= J Du Du^ 8{uiu a - 1) exp i J dr ^-u\ (zV M ) 



a 



Using the identity, 



we can present eq. fl2"5D in a surface form: 



_i ( „^„ ) + (Vill) t (V; 



(25) 



(26) 



W£^ N) = J DuDu ] 5{\u\ 2 - 1) expi J dS^ 



-{u ] F^ u u) + i{V^u) ] [V v u) 



(27) 



where is the field strength in the fundamental representation. Eq. (^) has been first pub- 



lished in ref.|T3| however with an unexpected overall coefficient 2 in the exponent. Eq. 
presents the non-Abelian Stokes theorem for the Wilson loop in the fundamental represen- 



tation of SU(N). In the particular case of the 577(2) group transition to eq. fll5D is achieved 
by identifying the unit 3- vector: n a = «j 1 (r a )^M /3 where 



1 2 



u 



cos I e" 
sm | e 2 



2iu^d T u = a (cos f3 — 1) + (a + 7). 
It should be mentioned that the quantity 



dadr i diU^djU 



2nQ 



(2t 



(29) 



appearing in eq. fl2"?D is the topological charge of the 2-dimensional CP N ~ 1 model. For closed 
or infinite surfaces Q is an integer. For example, in case of the 'Wu-Yang monopole' u a can 
be chosen in two forms compatible with periodicity: 




• §e*j ° r f°sin| )■ ^ 

They are regular everywhere except the negative (positive) z axis. One can regularize the 
string singularity of u a relaxing the condition \u\ 2 = 1. The coefficient 2ir guarantees that 
eq. ( p7|) does not actually depend on the choice of the surface even for singular it's, cf. the 
discussion around eq. ([H]) above. At the same time it means that CP"' 1 'instantons' of 
an effective 2-dimensional model, having integer topological charges, can hardly be relevant 
to the area behavior of the Wilson loop, as conjectured recently in ref.Q. The coefficient 
with which the topological charge enters the formula for the Wilson loop (B7|) corresponds 
to the 'instanton angle' 9 = 2ir, hence it is unobservable from the point of view of the 2- 
dimensional instantons. Only configurations with half-integer topological charges (like that 
given by eq. (|30D or its gauge equivalents) stand a chance of being of relevance to confinement. 

In case the Wilson loop is taken in the adjoint representation of the SU(N) gauge 
group the combination irtiHi in eq. (|J) is the highest root. Only group elements of the 
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form exp(ia>iHi) commute with this combination, belonging to the maximal torus subgroup 
[/(l) 7 ^ 1 . Hence, in case of the adjoint representation one in fact integrates over the maximal 
coset SU^/Uil)*- 1 = F N ~\ i.e. over flag variables |J, §. 



4 Reply to criticizm by Faber, Ivanov, Troitskaya and 
Zach 



In a recent paper |TTJ Faber, Ivanov, Troitskaya and Zach (FITZ) have questioned the va- 



lidity of our formula for the Wilson loop, eq. (f|). Their points can be summarized as follows: 

1. A direct calculation using eq. (^) in two simple cases where the Wilson loop is known, i.e. 
for a pure gauge potential and for a vortex field, produces zero results instead of correct ones. 

2. The regularized evolution operator for an axial top used by us to derive the formula for 
the Wilson loop, when computed directly, is zero. 

3. The evaluation of this evolution operator by introducing small regulator moments of in- 
ertia, as suggested by us, is prohibited because it changes the symmetry from SU(2) to that 
of[/(2). 

4. Another formula for the Wilson loop is proposed which is similar to but different from our. 

We consider these points below, one by one, and show that all four are groundless as due 
to errors in mathematics. 

1. In sections 4,5 of their paper FITZ attempt to compute the Wilson loop in two simple 
cases: i) in a pure gauge background and ii) in a background of a vortex. In both cases the 
Wilson loop is known. In this calculation FITZ use our discretization of the path integral 
though without the regulator terms which, as we have stressed in the original paper [|]], are 
important to get the correct result. Nevertheless, it could be a useful exercise, were it not 
for mistakes in mathematics. 

In this calculation FITZ use several relations for group characters citing a paper by 



Bars |T7| "modified for our case". The modification is not performed properly. Ref. 
deals with the GL(N) group, and averaging is performed over the U(N) group, not SU(N) 
which is the case in question. The characters of a group are unambiguously defined for 
the elements of that group. Meanwhile, FITZ intensively use such ambiguous quantities as 
SU(2) characters of the non- SU{2) elements like t 3 , t\ = 1/4, etc. (t 3 is one half of the 
Pauli matrix T3). Though no explicit definition of the characters of non-S77(2) elements is 
given in the paper, from eqs. (64) and (A3) of ref.[n one infers that FITZ implicitly use 



the definition: 

xAih) n U}= £ m^DUiU) (31) 
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where D^JJJ) are Wigner's finite rotation matrices, Em^4m(^) = Xit^l- [The properties 



of .D-functions are given, e.g., in ref. [18 



The key formula of FITZ calculation is their decomposition formula (51), 



exp zTr[t 3 U] = E Xj[h U] = J2 a oi z ) E mD Lm(U), 

j j m=-j 



(32) 



and its inverse (eq. (53) of fill), 



aj(z) 



dU Xj [hU ] ] exp(zTr[t 3 U]). 



(33) 



3(3 + 1) 

The quantities Xj\k?>U} do not form a complete set of functions for the space of SU{2) 
matrices, therefore there are no reasons why the decomposition (|32]) is at all possible, for 



whatever coefficients aj(z). To see that it is in fact wrong let us present the decomposition 



of the l.h.s. of eq. 
of the fact that £3 
element of SU(2). 
We have 



using well-defined characters of group elements only. We make use 
—i/2)(ir^) = (—i/2)exp(iirts) where the last factor is definitely an 



exp zTr[t 3 U] = exp {(-^/2)Tr[e i7r * 3 U}} = E 



a A z )Xj 



u 



E a," W E ^ D L m ( U )i aj(z) = e** (2j + 1) 



2J2j+i(z) 



m=-j 



(34) 



the coefficients aj(z) being well-known from the lattice strong-coupling expansion [^0 



Eq. (0) differs significantly from eq. ( |32"D suggested by FITZ; there exists no choice of 
coefficients aj(z) for which eq. (|32|) coincides with the correct one, eq. (0). Eq. (^) and its 
inverse, eq. (|33|), (eqs. (51) and (53) of [Tl[) result from uncritical 'modification' of ref.fTTH 
and are wrong. 

Furthermore, assuming the decomposition ( |3"2"D to hold true FITZ arrive to a "complete- 
ness condition" for their coefficients cij(z) (see their eq.(59)): 



+ Eoi(j + iH 2 (^ 

j>0 6 



(35) 



where (see eq.(A6) of JTT 



~Ji(z), 3 = 0, 

a ^ = \ ^1^' 3 = 1/2, 3/2, 5/2, 
0, 3 = 1,2,... 



(36) 



We plot the l.h.s. of (p5[) as function of z in Fig.l: it does not look as being identically 
unity, as claimed by FITZ. It is another manifestation of that eqs. (^,^) are wrong 

1 Ironically, FITZ have checked their "completeness condition" themselves but only at one value of z = 1 
where the l.h.s. of eq. ( p5| ) is still rather close to unity being equal to 0.986. Having computed this value 
summing up the series on the l.h.s. up to j — | FITZ note: "Thus, the series converges slowly to unity" not 
paying attention to that their J = § contribution is ' 
The series converges rather rapidly but not to unity. 



10 4 and the next contribution at j = I is ~ 10 
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Figure 1: The l.h.s. of the "completeness condition" claimed by FITZ to be identically unity. 

Unfortunately, this "completeness condition" together with the erroneous coefficients 
dj{z) are at the heart of the calculation of the Wilson loop by FITZ both for the pure gauge 
and vortex cases, which cannot be, thus, considered as correct. 

2. In section 6 of their paper FITZ attempt to compute the (regularized) evolution 
operator for the 'Wess-Zumino' action, following directly our approach. This calculation 
has been presented in some detail in the original paper [Q], however, FITZ seem to be 
dissatisfied by it and present their own. Their final answer (eqs. (98) and (112) of ref. |l ]]]), 
which differs from our, is a result of several mistakes. 

First, going from eq.(83) to eq.(87) FITZ use a strange relation, 



instead of the correct (and trivial) exp(i2(N + l)I±/5), where I± and 6 are constants and N 
is the number of pieces in which one divides the contour. 

Second, and more important, both equations in (91) are erroneous, they do not follow 
from eq.(89) from where they are derived. Passing from eq.(89) to eq.(91) one gets: 



FITZ have written these formulae without the crucial factors cos (3 n . Because of this mistake 
the subsequent integration over Euler angles a, 0, 7 becomes gaussian, and the evolution 
operator for the axial top, as computed by FITZ, in fact becomes that of a free particle, 
which is definitely wrong. The factors cos f3 n being reinstalled, the derivation returns to that 
of our paper [|l| . 




(37) 



Tr^i^+i) = 2 - ^ $ a l + fifin + $ll + 25a n 57 n cos/5 n , 

Tr(R n Ri +1 r 3 ) = i(6a n + 5j n cos(3 n ), 

Sa n = a n+1 - a n , 5(3 n = (3 n+1 - fi n , Sj n = 7 n+i - 7„. 



(38) 
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3. The last objection by FITZ is to our alternative (and in fact equivalent) derivation 
of the evolution operator, this time through the standard Feynman representation for the 
path integral as a sum over intermediate states. FITZ quote our result for the evolution 
operator of an axial top with the ' Wess-Zumino' term, evolving from its orientation given 
by a unitary matrix Ri at time t\ to orientation R 2 at time t 2 . 



Z^RiiR-y) = (2J + l)Djj(R 2 Rl) exp 



-i{t 2 -h) 



J 



2/, 



(39) 



where I± is a regulator moment of inertia, I± — *■ 0. Apart from a nontrivial dependence 
on the orientation matrices Ri 2 coming through the Wigner D-function, this expression 
contains a phase factor exp(— i(t 2 — ii)...). It is an overall factor independent of the external 
field: it can and should be absorbed into the integration measure to make the evolution 
operator unity for the trivial case R 2 = R\ = 1. Indeed, dividing the time interval into iV 
pieces of small length 5, N5 = t 2 — t±, one can write this factor as a product, 



exp 



J 

277 



N 

IJexp 

k=l 



45 



J 



2h 



(40) 



where, according to the regularization prescription of ref. 6/I± <C 1 so that each factor 
is close to unity Each factor can be now absorbed into the integration measure dR(tk) 
in the functional- integral representation for the evolution operator (|39|) . The fact that the 
factor is complex is irrelevant; moreover, it is typical for the path- integral representation of 



the evolution operators to have a complex measure, see the classical Feynman's book fl9 |. 
However, FITZ write: "...a removal of the fluctuating factor is prohibited since this leads to 
the change of the starting symmetry of the system from SU(2) to £7(2)". It may seem that 
FITZ believe that an absorption of a constant factor into the integration measure changes 
the number of variables over which one integrates. 



4. FITZ make an attempt to derive another formula for the Wilson loop (see section 2 
of ref. ||11|| ); this derivation is, however, inconsistent. Their resulting path- integral repre- 



sentation for the Wilson loop is (see eq.(30) of |L1]; we simplify their notations to make the 
formulae better readable): 



Wj = irnrWi I £ ( 2J + ^ exp V { drm ^ y/™[A n {T)A°(T)]) , (41) 

{ZJ -f- L) J {m(r)} V JC 7 

where A n = QAQ* 1 — iQQ^ is the gauge-transformed Yang-Mills field tangent to the con- 
tour parametrized by r, < r < 1, such that A = A^dx^ /dr, Q = d^VLdx^/dr. One 
integrates over all gauge transformations ^(t) i n a given representation of the SU(2) group, 
labelled by spin J, and sums over the projections m from — J to + J at all points of the loop. 

Integration over a variable m(r) assuming only integer or half- integer values is a unusual 
construction and one can question whether such an integral has a limiting value as one takes 
the number of discretization points iV to infinity. We have checked that eq. ( f41~D does not 
possess a finite limiting value at least for the zero field, A = 0. However, this is not our 
main point. The main problem is the derivation of eq. ( |4"1"D by FITZ. 
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If one looks into how eq. fl4"I|) has been derived by FITZ one finds that it is a result of 
two mistakes in mathematics. 

FITZ start from using the following discretized form of the Wilson loop (see their eqs. 
14,15): 



W r = — lim 

di n ^°° 



d rXr [fi(r 7V )(l + iAAr)fi t (r^_ 1 )] ...d rXr [n(ri)(l + %A At)&(t n )] , At = T k - r fc _i, 

(42) 

where gauge transformations Q and the gauge field A are taken in a given representation r 
with dimension <i r ; Xr is the group character. Indeed, integrating over gauge transformations 
at all points along the loop O(rfc) and using the relation 



dQ(Q\ lbl (Q) a2b2 



tL 



(43) 



one recovers the Wilson loop as a path-ordered product of the factors (1 + iAAt) along the 
loop. However at this point FITZ depart from a consistent derivation. 

It is tempting to say that, since the discretization points on the contour can be taken as 
close to one another as one wishes, the gauge transformations Q at neighboring points are also 
close. It should be emphasized that this is, generally, wrong: one needs a fully independent 
integration over fi's at neighbor points - otherwise the unity matrix on the r.h.s. of eq. fl4"3]) 
is not achieved. In other words the unitary matrix f2(rfc)f^(rfc_i) is, generally, not close to 
the unity matrix, even though the points Tk and r^-\ are close. The derivative Q does not 
exist in a strict sense since taking the neighbor points closer does not change the fact that 
^(77.) and fi(Tfc_i) are independent integration variables. 

This important circumstance is neglected by FITZ who say: "Due to the infinitesimality 
of the segments we can omit the path ordering operator" and write (see their eq. (18)) 



n(r k )(l +iAAr)^( 



T k -l 



exp 1 



drA i 



(44) 



This equation is erroneous for reasons explained above. Only if fi(r) has a finite derivative 
one can expand n(rjfe)0'(rfc_i) = 1 + £X T *^-i)^ ( r fc-i) At but then the exponent on the r.h.s. 
of eq. ( f44"D should be expanded too. If the derivative is large so that one cannot expand the 
exponent (and FITZ keep it), then eq. ( fPj) is simply wrong. To see it more clearly let us 



take an example of a zero field, A = 0, A 



iQQ^. Apparently, 



The path ordering on the r.h.s. 



(45) 

is absolutely 



if the l.h.s is an arbitrary unitary matrix 
necessary to restore the l.h.s. 

Even if one accepts the wrong eq. (HI) the way FITZ proceed further is questionable. 
According to eq. (|42|) FITZ need to compute the character of eq. (f44"l). Denoting the exponent 
in eq. (03) by an su(2) matrix $ in representation J, 



drti 



d:T 



A' 



(46) 
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the character can be written as 

j j 



Xj[expi$] = Yl exp imV® a ® a = Yl ex P im J J dTi[A n {r 1 )} a J dT 2 [A n (T 2 )] a 



m=—J m=—J 
.1 



£ exp imJ J j dnd^cjTriA^i^A^i^)], (47) 



m=—J 



J(J+1)(2J+1)' 

Meanwhile, FITZ use the following formula (their eqs.(19), (21) and (22)) 

Xj[expz$] = J2 ex P im J dr ^2Tr[A n (r)A n (T)}. (48) 

m=—J 

Leaving aside the incorrect numerical coefficient, the square root of a product of integrals is 
not equal to the integral of the square root of the product of integrands. In this way FITZ 
arrive to their formula for the Wilson loop, eq. (fllD . Their SU(3) generalization (section 3 
of |Tl|) is based on the same manipulation and is wrong from the start. 



To conclude this section: all points of the FITZ' criticism of our formula for the Wilson 
loop are based on their errors in mathematics, and the alternative formula suggested by 
these authors is mathematically inconsistent, too. 

Finally, we would like to comment on two sentences from FITZ' paper. In the Conclusion 
they write: "The use of the erroneous path integral representation for Wilson loops [meaning 
our formula] has led to the conclusion that at large distances the average value of the Wilson 
loops shows area-law falloff for any irreducible representation of SU(N)..." , and "...has led 
to result supporting the hypothesis of maximal Abelian projection". We are not aware of 
any work where either of the statements has been mathematically derived using our formula. 

Obviously on the contrary, the non- Abelian Stokes theorem of section 3 stresses the differ- 
ence between various representations for Wilson loops: in SU(2) there is a clear distinction 
between integer and half-integer representations (see, e.g., eq. (0)); in higher groups the 
number of integration variables itself depends on the representation (see the end of section 
3). 



5 Lattice-regularized formula for the Wilson loop 

The Wilson loop on a lattice is just the trace of the product of link variables in a given 
representation along the discretized contour. Else, it is the character of the product of link 
variables, 

Wj = 2J 1 + - X j (C/Ar,JV-l^JV-l,Ar-2 • • • U 1;N ) (49) 

(we consider the SU(2) gauge group for simplicity). The aim of this section is to derive a 
representation for eq. ( |49D analogous to the continuum eq. (f|), using lattice regularization. 
It has the form: 
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r N z 

Wj = AT 1 / J] dS k exp - Tr (^U^S^n), (50) 

k=l ~ 

where z is a function of J and M is a normalization coefficient, both of them to be determined 
below. Uk,k-i are link variables which in the continuum limit become Uk,k-i ~ ^+ia.A k where 
a is the lattice spacing and A k is the component of the Yang-Mills field along the link at 
point k. In eq. (|50"D one integrates over all gauge transformations S k at lattice sites k along 
the loop, with the condition that So = SV where N is the total number of links in the loop. 
Integration is over the Haar measure normalized to unity. We shall see that eq. ( |50|) is valid 
for large Wilson loop with number of links N ^> 1. Notice that on cubic lattices N is always 
even for closed loops. 

In the limit aA k <C 1 and of smoothly varying A k and S k eq. ( |50D coincides with eq. (Q) 
provided one takes z = 2 J. However, eq. (Q) has been derived from another regularization of 
the functional integral over gauge transformations, and there is no a priori reason to expect 
that in the lattice regularization z should be the same. 

Let us expand the exponent in eq. flBUD according to eq. fl34f ): 

Wj = N~ l j l[dS k J2a jk {z) E ^< mt (^M-i4-i), (51) 

fc=l h m k =-j k 

~ aj {z) = (2j + l)h 2j+1 (z). (52) 

Every matrix S k enters twice this expression. The integration can be performed as follows: 
dSk D 3 ^ imk+i (AS k ) D^ kmk {S\B) = — — - 5 jkjk+1 5 mkmk+1 D{^ kmk (AB) . (53) 

Using eq. (|5^) we get: 

Wj = Af-^MzT^-^EDLrniU^N-iU^ (54) 



J 



in 



2 

bj(z) = -J 2 j+i{z) (55) 

The factor e **W-™) = 1 for any j, m since j— m is an integer and N is even. Consequently, 

Wj = N' 1 fiMf Xj (U NtN - X UN-i^-2 ■ ■ ■ U 1>N ) , (56) 
j 

AT = E(2j + l)feWr, (57) 



where we have chosen the normalization factor M such that Wj = 1 for unity link variables. 

Thus, eq. ( |50"D is actually a weighted sum of Wilson loops in all representations j. However 
at N ^> 1 the sum in eq. (|56"D is dominated by the term with j maximizing the coefficient 
bj(z). Then only one term survives both in the numerator and the denominator of eq. 
and we obtain: 
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Wj 



iV>l 



2J + 



yXj {U N ,N -\U N -1,N -2 ■ ■ ■ U^n) , 



as it should be for the Wilson loop in representation J. 

Let us choose z(J) from the requirement that bj(z) 
for J < 2 one can choose 



is maximal at j 



(58) 

J. We find that 



*(J) = 2(J+1). (59) 

With this value of z(J) the Wilson loop in the needed representation J is reproduced by 
eq. (|50D at large number of links. [We would like to mention on this occasion that one can 
obtain eq. (|59"D also in a continuum formulation where the regularization is performed by 
introducing small moments of inertia I_ui, as in ref. 0. In that paper we used the following 
limiting procedure: first 7ii — > 0, then I± — > 0, and obtained z(J) = 2J. Had we chosen the 
opposite order of limits we would get eq. (|59T).1 

In fact, one is not bound to take simple values like given by eq. (|59| ) but choose the values 
of z(J) for J = |, 1, | . . . such that bj(z(J)) is maximum maximorum at j = J. The result of 
such 'fine tuning' is presented in Table 1, where we give the best value of z for given J, the 
maximum value of &best and the ratio of bj(z) to &best for the most 'dangerous' competitor at 
./ •/• 

Table 1 



J 


best z 


&best 


bj(z)/b hest (j) 


1 

2 


3.25103 


0.296 


0.50 (0) 


1 


4.36765 


0.198 


0.60 (1/2) 


3 
2 


5.46564 


0.146 


0.86 (0) 


2 


6.55104 


0.114 


0.83 (1/2) 


5 

2 


7.62728 


0.093 


0.83 (3) 


3 


8.69644 


0.078 


0.85 (7/2) 



If, for example, we wish to get the Wilson loop in representation J = |, we choose the 
coefficient z(5/2) = 7.62728. The largest contribution to eq. (|56|) will be from j — J — |. 
The next largest but parasite contribution will be from j = 3; however, the Wilson loop in 
this representation will be suppressed as 0.83^ as compared to the needed j = | contribution. 
If the number of links in the loop is iV = 16 the suppression factor is 0.051, for TV = 24 it is 
0.011. 

The lattice version of the 'non-Abelian Stokes theorem' is trivial: one takes any surface 
(realized as a collection of plaquettes) inside a given discretized contour (realized as a chain 



Sk Uk,k- 



Sk-l ^3 



for each link 



of links) and writes in the exponent a sum of the terms Tr 
belonging to all plaquettes inside the contour. Since all internal links are encountered twice 
in this sum, once going in a positive and once in a negative direction, all links cancel except 
those lying at the border of the surface, that is on the contour itself. The cancellation of 
links is due to the fact that Tr(\^rs) = — TriVWs) for a unitary matrix V. The lattice version 
of the non-Abelian Stokes theorem reads: 
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Wj = M- 1 Jll[dS k \e X pl~ E Tr{StU k)k ^S k ^T 3 ) . (60) 

yfcGsurf / plaqSsurf linksGplaq J 

We would like to stress that any Stokes theorem is trivial when written in discretized form. 
At the same time, if we first assemble links belonging to one plaquette and ascribe to the 
quantities S and U arguments corresponding to the plaquette centers, we recover, in the 
appropriate continuum limit, the non-Abelian Stokes theorem for the continuum, eq. (|T5|). 

If one uses eq. (|50|) or eq. (^) for the average Wilson loop in the lattice formulation of 
gauge theory there is no need to perform explicitly integration over gauge transformations as 
these are included in the integration over link variables. Therefore, one can use the following 
formula for the average Wilson loop on the lattice: 



(Wj) = M 1 \\ I dUi exp [lattice action] exp 
i J 



9 

^ lec 



(61) 



Actually, this formula has been suggested in ref. ||, however with the coefficient z(J) = 2 J 
arising from another regularization. This is not too consistent: in lattice calculations one 
has to use lattice regularization for the Wilson loop. Therefore, one has either to put 
z (J) = 2( J + 1) (for small J) or, better, pick it up from Table 1. 

Finally, we note that the choice of the matrix T3 in eq. ( |BTD is arbitrary: one can take 
any rotated matrix as well. This fact can be used to increase many times the statistics in 
numerical computation of the average Wilson loop from eq. (|61|). It may be interesting to 
study the dependence of the r.h.s of eq. (|61~D on the coefficient z for a fixed but large loop as 
it can provide new information on the difference of integer and half-integer representations 
in respect to confinement. 



6 Conclusions 

We have formulated the non-Abelian Stokes theorem for the Wilson loop. The path-ordering 
is replaced by an ordinary exponent of a surface integral, but at the price of an additional 
functional integration over all gauge transformations of the Yang-Mills potential, actually, 
over a coset depending on the representation in which the Wilson loop is considered. We have 
given several forms of this representation, and discussed requirements on the continuation 
of the fields inside the contour. 

Since the validity of our formula for the Wilson loop, which is key to the non-Abelian 
Stokes theorem, has been questioned recently in ref. [|ll[] we had to reply to that criticism. 
We have demonstrated that it is groundless as due to mistakes in mathematics, which we 
have thoroughly pinpointed, one by one. 

As the lattice regularization of functional integrals is one of the most popular we have 
included the derivation of the formula for the Wilson loop in lattice regularization. The 
resulting eq. (|50| ) is very similar to the continuum formula (^). The corresponding lattice 
version of the Stokes theorem is almost trivial. 

A formula analogous to the non-Abelian Stokes theorem can be also derived in gravity 



theory for holonomies in curved spaces [21 
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